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ABSTRACT 

We present a Chandra X-ray observation that shows, unambiguously for the first time, the presence 
of a giant diffuse X-ray-emitting corona around the edge-on disk galaxy NGC 4631. This corona, with a 
characteristic thermal temperature of 2-7 x 10 6 K, extends as far as 8 kpc away from the galactic plane. 
The X-ray morphology resembles the radio halo of the galaxy, indicating a close connection between 
outflows of hot gas, cosmic rays, and magnetic field from the galactic disk. Enhanced diffuse X-ray 
emission is apparently enclosed by numerous Ha-emitting loops blistered out from the central disk of 
the galaxy, as is evident in a comparison with our deep Hubble Space Telescope imaging. 

Subject headings: galaxies: individual (NGC 4631) — galaxies: ISM — galaxies — spiral — X-rays: 
galaxies 
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1. INTRODUCTION 

Though postulated more than four decades ago (Spitzer 
1956), the existence of a hot (T ~ 10 6 K) extended 
corona around our Milky Way galaxy is still subject to de- 
bate. What limits our ability to understanding the corona 
around the Milky Way is our location in its midst. Given 
the expected characteristic temperature of ~ 10 6 degrees 
for this gas, soft X-ray emission is the primary tracer for 
the Galactic corona. However, the soft X-ray sky observed 
from our location within the Galaxy is very confusing due 
to the multiple origins of X-ray-emitting hot gas. Any ex- 
tended Galactic corona, if present, is embedded between a 
nonuniform foreground, produced by the Local Hot Bub- 
ble around the Sun and various contributions from other 
nearby interstellar gas clouds, and an extragalactic soft X- 
ray background of an uncertain spectrum. Disentangling 
these components along a typical line of sight is extremely 
difficult, if not impossible. 

The most direct way to study extended galactic coronae 
is to observe nearby edge-on disk galaxies. ROSAT ob- 
servations of such galaxies have demonstrated that extra- 
planar X-ray emission is present in relatively normal spi- 
rals (NGC 891, Bregman & Houck 1997; NGC 4631, Wang 
et al. 1995). Such "normal" edge-on galaxies, in which star 
formation is not dominated by galactic nuclear regions, but 
extend over a large fraction of the galactic disk, are ideal 
for investigating general properties of galactic halos and 
their connections to galactic disks. The limited sensitivity 
and spatial resolution of ROSAT, however, did not allow 
for an unambiguous separation of diffuse emission from 
discrete sources, let alone a detailed study of morphology, 
ionization, and dynamics of hot gas. 

In this Letter we report preliminary results from ob- 
servations of the SBc/d type galaxy NGC 4631 with the 
Chandra X-ray observatory and the Hubble Space Tele- 
scope. This galaxy is ideal for studying soft X-ray emission 



from its galactic corona because of the very low Galac- 
tic foreground absorption (iVni = 1.2 x 10 20 cm -2 ) and 
its edge-on orientation with respect to the viewing direc- 
tion. Active star formation is apparent throughout much 
of the NGC 4631 disk, apparently triggered by interactions 
with neighboring galaxies. Evidence for this interaction in- 
cludes the presence of HI tidal tails (Weliachew, Sancisi, 
& Guelin 1978). Such active star formation is presumably 
a necessary condition for the heating of coronal gas. 

2. CHANDRA OBSERVATION AND DATA CALIBRATION 

Our Chandra observation was carried out with the Ad- 
vanced CCD Imaging Spectrometer chip S3 (ACIS-S3) at 
the focal plane of the telescope mirror. The spatial res- 
olution of the observation before any smoothing is ~ 1", 
corresponding to ~ 37 pc at the 7.6 Mpc distance of the 
galaxy. This superb resolution, together with the high sen- 
sitivity of the instrument and the long exposure (~ 60 ks), 
allows us to remove point-like sources and to directly probe 
the corona and its relation to other galactic components. 
We calibrated the data using primarily the software pack- 
age CIAO v2.1 with the latest calibration products (caldb 
version 2.3). We selected data only in the energy range 
of 0.3-7 keV. Particle background is greatly enhanced at 
higher energies, while the calibration of both the gain (un- 
certainties up to ~ 30%) and quantum efficiency (up to 
~ 10%) of the detector are unreliable at lower energies. 
We excluded bad pixels, rows and node boundaries, and 
screened the data for high and low background times by 
clipping observation periods with an offset exceeding ±3<r 
from the mean quiescent count rate (1.16 cts s~ x ). This 
filtering led to the exclusion of 8% of the total counts and 
6% of the total exposure time, resulting in a cleaned ex- 
posure time of 55.1 ks. 

We further conducted both background subtraction and 
exposure correction of the data. We use deep (138 ks) 
blank-field background datasets with similar foreground 
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absorption and identical focal plane temperature as our 
NGC 4631 observation to determine the position depended 
background rates. The same data filtering criteria were ap- 
plied. We created bad pixel and node-boundary removed 
exposure maps in various energy bands. For the spec- 
tral analysis of extended X-ray emission, we constructed 
weighted instrument effective area and spectral response 
matrices, using a low-resolution (bin size = 4"), 0.3- 
2.5 keV band image of diffuse emission in detector co- 
ordinates. The weighting corrects for variations (up to 
~ 20%) in the instrument sensitivity across the detector 
(e.g., effective area, gain, and quantum efficiency). It is 
important to note that the currently available instrument 
response matrices at low energies (< 1 keV) employ var- 
ious approximations, which may not be appropriate and 
may lead to underestimation of line strength. Since NGC 
4631 fills nearly the entire ACIS-S3 field-of-view, we ex- 
tracted the background from the same detector regions in 
the empty-field background observations described above 
as the source regions in the NGC 4631 observations. 

3. ANALYSIS AND RESULTS 

Chandra images in the soft and hard X-ray bands are 
presented in Fig. 1. Whereas there is little indication 
for extra-planar emission in the hard band, the enhanced 
emission in the soft band extends as far as 8 kpc away from 
the plane of the galaxy, preferentially towards the north. 
The diffuse emission becomes progressively more extended 
with decreasing energy (Fig. 2). The distribution of the 
low surface brightness emission is nearly symmetric, rela- 
tive to the minor axis of the galaxy. This extra-planar soft 
X-ray emission clearly represents the corona of the galaxy. 

We have conducted a preliminary X-ray spectral anal- 
ysis of the corona. We extracted a point-source excised 
spectrum in a corona region 4' wide and 2'5 high (cf. Fig. 
2b). The spectrum in the 0.3-2 keV range can be char- 
acterized reasonably well (x 2 /d.o.f. = 78/47; Fig. 3) by a 
thermal plasma model with two temperature components: 
0.18(0.16-0.20) keV and 0.61(0.48-0.71) keV. The jointly 
fitted metal abundance and absorption column density are 
8(5-12)% solar and 3.7(2.2-5.8) x 10 20 cm" 2 , which is con- 
sistent with the expected foreground absorption. The un- 
certainties quoted in the parentheses are 90% confidence 
intervals. Assuming that the model is a reasonable char- 
acterization of the corona emission, we derive a conversion 
ratio of the emission measure (in units of cm~ 6 pc) to 
the observed count intensity (10 -4 counts s _1 arcmin~ 2 ) 
of 91 (27) for the low (high) temperature component. The 
average ratio of the low to high temperature emission mea- 
sures is ~ 8 over the corona region. The low-temperature 
component dominates the outer corona (> 2'), whereas 
the high-temperature component becomes significant in 
regions close to the disk. The characteristic gas density 
and thermal pressure (p/ks) are ~ 1(5) x 10~ 3 cm~ 3 and 
~ 5(40) x 10 3 cm~ 3 K in the outer (inner) corona. The 
total luminosity of the corona is a few times 10 39 ergs s _1 
in the 0.3-2 keV band; the uncertainty in the absorption 
correction is large. More careful analysis is required, espe- 
cially in regions close the disk. Nevertheless, the amount 
of observed diffuse X-ray emission clearly accounts for only 
a small fraction (< 1%) of the estimated energy input from 
supernovae in the galaxy (Wang et al. 1995). 



The energy-dependent distribution of the diffuse X-ray 
emission is further illustrated in Fig. 4. The emission in 
the lowest energy range (0.3-0.6 keV) has the flattest dis- 
tribution. The dip just below the galaxy's major axis 
(z = 0'), most prominent in lower energy bands, is caused 
by X-ray absorption from the galactic disk which is slightly 
tilted with its near side below the major axis. The dip at 
z = 2' spatially coincides with a giant dusty arch discov- 
ered recently by Neiningcr & Dumke (1999). In fact, this 
arch seems to enclose the enhanced soft diffuse X-ray emis- 
sion, approximately outlined by the outer white contour in 
Fig. 2a, in the northern central region of the galaxy. 

4. MULTIWAVELENGTH COMPARISON 

Fig. 2b presents a morphological comparison of the dif- 
fuse X-ray and Ha emissions from NGC 4631. Whereas 
the Ha line emission represents primarily warm gas photo- 
ionized by stars, the soft diffuse X-ray emission arises in 
coronal gas. The coronal gas is distributed further out 
from the disk than the warm gas, although the galaxy 
might also have a very low surface brightness Ha-emitting 
halo, which may extend further out (Donahue, Aldering 
& Stocke 1995). The enhanced X-ray emission arises in 
regions containing coherent vertical Ha filaments, primar- 
ily in the upper central portion of the galactic disk. The 
left one quarter of the disk, though highly disturbed, does 
not show many vertical Ha filaments or much enhanced 
extra-planar diffuse X-ray emission. 

This association between enhanced Ha/X-ray filamen- 
tary features appears most vividly in Fig. 5, which covers 
the inner 5.5 kpc x 5.5 kpc region of the galaxy (cf. Fig. 
2b). We obtained this deep Ha image using the Wide- 
Field Planetary Camera 2 (WFPC2) on board the Hubble 
Space Telescope (HST): eight orbits narrow-band imaging 
plus one broad-band continuum. The image shows numer- 
ous filaments emanating from the disk, mostly apparent in 
the upper-right part of the image. Such coherent loop-like 
features most likely represent "a froth of merged superbub- 
bles", which were created by massive star clusters, have 
broken out from the dense neutral gas of the disk, and are 
expanding into the halo (e.g., MacLow, McCray & Norman 
1989; Norman & Ikcuchi 1989). The overall correlation be- 
tween this froth and the enhanced diffuse X-ray emission 
is apparent. However, relatively bright loops do not seem 
to have a configuration with open tops, as one might ex- 
pect in galactic chimney models. The prominent V-shaped 
double filaments, just above the large dark cloud near the 
center of the image, may not be an opened chimney and 
most likely represent walls of many projected loops. X-ray 
emission within the loops is generally enhanced by a factor 
of up to ~ 3. Furthermore, the WFPC2 image reveals sev- 
eral bow-shaped filaments far away from the central plane 
of the galaxy. The largest and brightest one, as labeled 
as "Arc" in Fig. 5, "covers" the top of the V-shaped fila- 
ments. These Ha filaments may represent the illuminated 
edges of in-falling clouds, which could originate in cooled 
halo gas or in instabilities of the dusty arch. 

The overall X-ray morphology of the corona (e.g., Fig. 
lb) resembles the well-known radio halo of this galaxy 
(Hummel & Dettmar 1990), indicating a close link between 
outflows of hot gas and cosmic ray/magnetic field from the 
galactic disk. One scenario for the formation of such a ra- 
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dio halo is that hot gas draws magnetic fields (and cosmic 
rays) from the disk into the halo, before being ultimately 
confined by the magnetic field tension at large off-plane 
distances (Wang et al. 1995). 

The unusually low metal abundance inferred from our 
spectral fit is apparently an artifact, which could be caused 
by the poorly calibrated instrument response shape be- 
low <~ 1 keV (§2) and/or by the overly-simplified spec- 
tral model we have assumed. A realistic temperature 
distribution should span a range, which tends to smooth 
out temperature-sensitive spectral features resulted from 
metal emission lines. While our spectral analysis provides 
an estimate of the temperature range, a realistic model is 
yet to be developed for a reliable measurement of the abun- 
dance. More importantly, however, the out-flowing gas 
may not be in an ionization equilibrium state (Breitschw- 
erdt & Schmutzlcr 1999). We are currently conducting 
a spatially-resolved spectroscopic analysis of the Chandra 
data and are investigating a self-consistent dynamical and 



thermal plasma emission model. 

Now with the corona around NGC 4631 firmly detected, 
one may naturally ask questions such as: does the extent 
of galactic coronae always coincide with the radio halos? 
Is the galaxy-galaxy interaction necessary to generate the 
coronae? Do properties of diffuse hot gas depend on the 
Hubble type of a galaxy? Ongoing and future observations 
with the Chandra and XMM-Newton X-ray observatories 
will allow us to answer these questions, leading to a better 
understanding of our own Galaxy as well as the structure 
and evolution of galaxies in general. 

We thank M. Markevitch and A. Vikhlinin for helping 
with the calibration software and the referee for useful 
comments, especially on existing problems with the spec- 
tral response matrices. This work is funded by SAO Chan- 
dra grant GO0-1150X, STScI grant GOO-8166.01, and 
NASA LTSA grant NAG5-8999. 
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Fig. 2. — Diffuse soft X-ray and Ho intensity distributions of NGC 4631. The Chandra image (a) is produced in the 0.3-0.9 kcV band 
and discrete sources (marked by crosses) are excised. The ground-based Ha image (b) is reproduced from Hoopes et al. (1999). The X-ray 
contours are at 1.3, 2.5, 4.2, 6.6, 10, 15, 21, 31, 45 X 10~ 4 counts s" 1 arcmin . The rectangular box marks the region from which the corona 
spectral data in Fig. 3 was extracted, whereas the square box illustrates the location of the inner region shown in Fig. 5. 

Fig. 5.— HST WFPC2 Ho image of the central 5.5 kpc X 5.5 kpc region of NGC 4631 (Fig. 2b). High level X-ray contours in Fig. 2 are 
plotted and a few tentatively identified Ha-emitting loops are outlined in right-hand panel. 
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Fig. 1. — Chandra ACIS-S images of NGC 4631 in the 1.5-7 keV band (a) and the 0.3-1.5 keV band (b). These images, both background 
subtracted and exposure corrected, are adaptively smoothed with the program CSMOOTH implemented in CIAO; the signal-to-noise ratio of the 
smoothed image is ~ 3. The far-UV intensity contours outline the morphology of the galactic disk. The field-of-view is 18 kpc X 18 kpc. 
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Fig. 3. — Chandra ACIS-S spectrum of the NGC 4631 corona. The histogram represents the best-fit two-temperature thermal plasma 
model (see text). 
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Fig. 4. — Diffuse X-ray intensity distributions along the minor axis of the galaxy NGC 4631. The profiles were calculated in a central 2' 
cut perpendicular to the galaxy's major axis (86° eastward of north). The contribution from discrete X-ray sources marked in Fig. 2a was 
excised. The northern side of the galaxy is to the right. 
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